Motion-reversal in a simple microscopic swimmer 
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We study the motion of a microscopic swimmer composed of a semiflexible polymer anchored at the 
surface of a magnetic sphere using hydrodynamic simulations and scaling arguments. The swimmer 
is driven by a rotating magnetic field, and displays forward and backward motion depending on the 
value of the rotational frequency. In particular, the system exhibits forward thrust for frequencies 
below a critical frequency uj* , while above uj* the motion is reversed. 

PACS numbers: 47.63.mf, 45.40.Ln, 87.19.St, 87.17.Jj 



Since the discovery of the microscopic world, scientists 
have tried to understand the directed motion of the or- 
ganisms that live at this scale. In particular, it has been 
found that many of the schemes for propulsion at low 
Reynolds numbers rely on beating or rotating flexible fll- 
aments periodically [TJ[5. In the case of rotation, the fll- 
aments typically have helical shapes (e.g. flagellae), and 
move through the fluid like a cork-screw in the direction 
dictated by the sense of rotation. Further studies in this 
area showed theoretically [3] |4] and experimentlally 
that simply rotating a straight semiflexible filament at 
its base is enough to achieve rectified motion along the 
axis of rotation. Within this propulsion scheme, the di- 
rection of motion turns out to be always from tail to head 
(i.e. the base). On the other hand, it has been recently 
shown by Dreyfus et.al. [6] that an actuated magnetic 
filament attached to a cell moves in the opposite direc- 
tion, i.e. from head to tail. These results suggest that 
there is still much to be learned about the coupling of 
the hydrodynamic and elastic interactions in model mi- 
croscopic swimmers, to which we refer more simply as 
microswimmers . 

Here, we introduce a novel microswimmer that moves 
in both directions depending solely on the value of the 
driving frequency. Our swimmer is composed of a semi- 
flexible filament (tail) that is anchored at the surface of 
a spherical magnetic bead (head). The swimmer is then 
driven by a rotatory magnetic field (see Fig. [T]). We 
have chosen this design because it is simple and could 
in principle be assembled in the laboratory. For exam- 
ple, one could self-assemble this system by using Strepta- 
vidin coated magnetic beads together with Taxol stabi- 
lized Biotinynlated microtubules [7], or it can be created 
at the lab scale using a similar approach as that em- 
ployed by Hosoi and co-workers to study low Reynolds 
number propulsion [5] . Also, our design is versatile and 
can be eventually used for pumping[5], as well as for local 
mixing in microfluidic chambers, and does not have the 
drawback of using live non-controllable microorganisms 

The origin of the forward motion displayed by our 
model microswimmer is due to the force that the fila- 



ment produces at the anchoring point, rectified by the 
magnetic torque acting on the sphere. The former force 
is a result of the bending stiffness of the polymer, which 
maintains a stationary deformed shape in order to min- 
imize the drag from the fluid. The mechanism responsi- 
ble for the switching dynamics from forward to backward 
motion is an instability that occurs when the magnetic 
moment lags the field vector by a phase (/> > 7r/2. This 
occurs when the driving frequency is sufficiently large 
such that the maximal torque applied to the system is 
not enough to overcome the drag on the rotating sphere 
+ filament. At this frequency, the magnetic dipole starts 
displaying large periodic oscillations along the direction 
perpendicular to the plane of rotation that give rise to a 
mix pattern of beating and rotation of the filament such 
that overall the system exhibits reverse motion. 

We study the motion of the proposed microswimmer 
by means of hydrodynamic simulations on a discretized 
version of our model. In particular, the magnetic bead 
of radius R and permanent magnetic dipole moment M 
is assembled using a set of 100 spherical particles of ra- 
dius OslIO]. The semiflexible filament of length L and 




FIG. 1: The proposed microswimmer is composed of a mag- 
netic bead (red sphere) with a semiflexible polymer (blue 
cylinder) anchored at its surface. The bead carries a per- 
manent magnetic moment M (depicted as the green arrow), 
and is driven by an external magnetic field B (gray arrow). 
Inset: Look of the swimmer after the discretization of the 
different components. 
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persistence length Ip is also discretized into N different 
particles of radius a / (see inset of Fig. [T| . For simplicity, 
we assume ag = af = a, which does not affect the overall 
behavior. The polymer is anchored by one of its ends to 
the south (or north) pole of the sphere, and held in place 
by harmonic springs on the first two beads. We have con- 
sidered this set-up in order to simplify the large design 
space, yet we note that other anchoring points or orien- 
tations may well give rise to an even richer dynamical 
landscape. 

The dynamics of the system is governed by 



(^) 0.75 



dt 



r, = ^/x,^..F,(t), 



(1) 



where is the position of the zth bead, ^l^ is the 
mobility matrix between all the particles in the sys- 
tem, and is the total force acting on the jth bead 
at time t. The mobility matrix /x^^ is taken to be 
the Rotne-Prager tensor. The force on each of the 
particles composing the sphere consists of two terms: 
F = 'Fu + 'Fb- The former corresponds to a harmonic 
force acting on each small bead, and can be written as 
F(7 = — Vrt/s, where Us = kJ2{fij — dij)^. In the last 
expression, r^j denotes the distance between particles 
i and j, and dij corresponds to the "equilibrium" dis- 
tance between both of these particles [TU]. The spring 
constant is set to fc = 2QQkBT / a? . The magnetic force 
Fb = 7B X M X r J /i? arises from the torque on the sphere 
due to the magnetic field that we consider to be given by 
B = J5(0, sin(ci;t), cos(wt)), where B is the magnitude of 
the field, and uj is the driving frequency. The numerical 
prefactor 7 is chosen such that the overall torque on the 
sphere is equal to B x M. In the case of the filament, the 
force on each of the particles results from changes in the 
potential energy described as Ufu = Us + Ub, where Us 
is a harmonic potential between adjacent beads with an 
equilibrium distance of 2a, and a spring constant equal 
to that of the springs linking the particles on the sphere. 
The bending potential is Ub — ~~ cos6'i), where e 

measures the stiffness of the filament, and 6i is the an- 
gle between adjacent bond vectors — r^+i — r,. The 
persistence length follows as Ip = ae/ksT. Finally, we 
re-scale all the lengths by R, the forces by ksT/a, and 
the time by the characteristic diffusion time of the sphere 
Tfi = B? /{iiRksT), where (Xr = l/(67r77i?) is the Stokes 
mobility of the magnetic bead. The dimensionless param- 
eters that govern the dynamics of the microswimmer are 
then the rescaled length of the filament L = L/R, the di- 
mensionless persistence length Ip — Ip/ R, the maximum 
magnetic torque (in units of ksT) otmax = BM/kBT, 
and the dimensionless driving frequency w = ujtr. 

In Fig. [2] (a) we show five different trajectories in the 
X direction for a set-up characterized by amax ~ 9500, 
L w 8.3, and Ip w 160. As noted before, the swimmers 
propels forward (positive slope) for values of u; < lu* , 




(b). 



0.02 0.04 0.06 0.08 



1 

0.5 
^ 
S -0.5 

q-H -1 
O 

^ 1 

•s 

^ 

0.5 

-0.5 
-1 



,\ .Ml ,A 'h 'A /A /A :l\ 'I') «=io8o 



0.004 0.006 



'T^ ' 
li \ Ml I i \ ' » \ I 



i\\ii\\ih\'J\'J\\i h\' 



0.002 0.004 0.006 0.008 



0.01 



u 






W4j 


r ] 




i\ to =2 


160 










/y\ yrv 









0.002 0.004 0.006 0.008 



0.01 



FIG. 2: (a) The distance (in units of the sphere radius R) 
traveled by the microswimmer in the x direction as a function 
of the rescaled time t/rii. The magnetic field and moment 
are held fixed such that the maximum torque a,nax ~ 9500. 
The filament length is L « 8.3, and the persistence length 
is given by Ip « 160. The different solid curves correspond 
to li = 1080,108,54,1620,2160 (from top to bottom). The 
dot-dashed curves correspond to the diffusive motion of the 
magnetic bead. Insets: amplification of the trajectories for 
the cases a) = 1080(upper inset), and 2160 (lower inset), (b) 
Time sequences of two components of the permanent dipole 
moment M of the sphere at three driving frequencies (indi- 
cated in the upper right corner of each graph). The thick solid 
curves correspond to the x component Mx and the dashed 
lines to the y component My. The thin solid line denotes 
the value of By{t)/B. All the sequences shown here have the 
same parameters as in part (a). 

which for this particular case corresponds to lu* ~ 1350 
(see below). For frequencies above w* the velocity di- 
minishes considerably and starts displaying reverse mo- 
tion (negative slope). Notice that the trajectories are 
smooth for uj < oj* , but above this frequency the system 
displays periodic forward and backward motion, that re- 
sults in a net negative speed if w is high enough (see 
lower inset). The origin of this behavior is due to the 
fact that the sphere not only rotates, but also "rocks" 
on the filament, and this produces a beat-like pattern. 
Actually, we define the value of w* as the point at which 
the swimmer starts displaying beating motion, which is 
directly related to the appearance of " strokes" of the per- 
manent magnetic dipole along the x axis. To illustrate 
this events, we present time sequences of M at frequen- 
cies below, at, and above the transition in Fig. [2] (b). 
From these sequences we first notice that the y compo- 
nent My of the dipole moment follows the driving field 
up to the transition point (upper and middle sequence). 
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At this particular value of w = , the x component of 
the permanent magnetic moment Mx starts displaying 
large strokes at intervals quite larger than the period of 
oscillation of the field (see middle sequence). This be- 
havior is retained for frequencies larger than lo* , except 
that the strokes become more frequent (lower sequence). 
Also, note that My seems to be a periodic function with 
a period larger than 1/Cb modulated by the field. The 
structure of My is nevertheless complicated and cannot 
be written in simple terms. However, the onset of the 
complex dynamics is found to correlate well with the ap- 
pearance of large amplitude periodic strokes in M^ ■ 

To rationalize this behavior, and in particular estimate 
the critical frequency lo* , we introduce a simple model for 
the microswimmer (see Fig. |3](a)). The most important 
assumption within our rationalization is that the dynam- 
ics of the system is stationary, and thus the deformation 
of the filament shape is constant in time. By minimizing 
an "energy" functional T = £ + W that corresponds to 
the sum of the deformation energy £ plus the work done 
on the system per unit cycle W we calculate the deforma- 
tion length I . Finally, by comparing the maximum torque 
with the dissipation of the sphere-ffilament system we 
obtain the critical frequency lu* . The energy of defor- 
mation within our model is given by £ ~ IpR^ksT /l^ , 
where we have considered that the deformed piece of 
the filament can be described by two arcs of a circle 
of radius r each spanning an angle 9. This simple re- 
sult captures to leading order in I the deformation en- 
ergy. The work done on the filament per cycle is given 
(to leading order in /) by W ~ u!T±kBTR^l/a'^, where 
T± = £^±a^ /ksT, and is the perpendicular drag con- 
stant of a filament segment of length a. Minimizing 
yields the deformation length I ^ a{lp/{aujTj_)y/'^ . The 
critical frequency is then simply calculated by setting the 
maximum torque BM equal to the dissipation of the ro- 
tating sphere+filament (,jfUj + W. This results in the 
following equation for lo* 

* / * \ 3/4 

where lo'^ — BM /^^ is the critical frequency for a sin- 
gle magnetic bead with no grafted filament, and is 
the rotational drag coefficient of the sphere. Note that 
this equation depends only on a single dimensionless pa- 
rameter Sn ~ R/a{Lo'^T±a/lp)^^^ that corresponds to a 
characteristic " sperm number" . This parameter can be 
conveniently written also as Sr ~ {amax/lpY^^ ■ The be- 
havior of LO* /loI as a function of Sr is shown in Fig. [s] 
(b). As expected, for large values of Sr (low rigidity or 
large magnetic fields) lo* /lo'}^ is essentially 1, while for 
values of Sr < 1 the critical frequency starts deviating 
considerably from lo'^ and can be easily shown to scale 

as LO* ^ ^*rS'^jI^ ^ Lo*j^{amax/lpY^'^ ■ Nevertheless, it 
should be noted that the dependence on the physical pa- 
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FIG. 3: (a) Model of the stationary state of the microswim- 
mer. The shape of the deformed part of the filament is de- 
scribed by two pieces of a circle of radius r. The overall 
deformation length is and 7? is the radius of the sphere, 
(b) The dependence of the rescaled critical frequency lo* jij*^ 
as a function of the characteristic sperm number 5*^ ~ 
Rjaiufj^T^ajly)^^'^ . The critical frequency in the absence of 
the filament is defined as ojj = EMj^r- 

rameters is very weak. Also, it is important to point out 
that we have neglected all the numerical prefactors in 
our derivation since our theory only depends on a single 
parameter whose functional form we obtained. The pref- 
actor for Sr, will in general be system dependent, and can 
be calculated by computing lo* /lo*j:^ from the simulations 
(or experiments) at a fixed value of OLmax, as is done in 
what follows. 

Finally, in Fig. |4] (a) we present the variance of the 
X component of the magnetic dipole a — {{M^ / M^) — 
{{Mx/M))'^ as a function of lo/lo*. The variance a is 
very sensitive to the appearance of oscillations in M^, 
and could be measured experimentally using induction 
techniques. The critical frequency lo* for the different 
parameters employed here is calculated from Eq. [2] using 
Sr — 4:-S{amax/lpy^^ ■ The prefactor is obtained from 
the data at amax ~ 9500, L « 8.3, and Ip « 160. It 
can be clearly observed from this graph that the mag- 
netic moment suddenly starts "oscillating" in the x di- 
rection with a large amplitude since a jumps sharply at 
the critical frequency. For larger frequencies, the ampli- 
tude of these oscillations becomes smaller, and hence the 
decrease in a. Interestingly, all the curves collapse when 
plotted in this way, implying that the dynamics of the os- 
cillations only depends on lo/lo^. In Fig. [4](b) we plot the 
corresponding velocites in the x direction v^/vd for the 
same parameters as in part (a). The velocity Vd — R/tr 
corresponds to the diffusive speed of the magnetic bead. 
As can be observed from this graph, the velocity grows 
as a function of the frequency for lo < oo* . Near (but 
below) lo* it essentially plateaus, and finally above the 
critical frequency it undergoes an abrupt decline in its 
value until negative velocities are attained. Interestingly, 
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FIG. 4: (a) The variance cr = {{M^/Mf) - {{M^/M)f, and 
(b) the rescaled velocity Vx/vd (see text for the definition 
of Vd) as a function of the rescaled driving frequency u/uj* . 
The different parameters used in each curve are displayed in 
the legend. The horizontal dashed line corresponds to the 
maximal achieved velocity w4™-e ~ 16wd- The vertical dashed 
line corresponds to the critical frequency at which the system 
becomes unstable (see text for details). 

as the frequency is tuned to even higher values the swim- 
mer starts moving forward again. The maximum for- 
ward velocity found for the parameters employed here is 
^iriaf ~ 16wd, for the case ctmax ~ 9500, L = 8.3, and 
Ip — 800. In a real set-up, these parameters would be 
equivalent to having a microtubule of length L ~ 10/i 
grafted on a 1/i radius bead, and swimming at a speed of 
ti ^ 8/i/s. In the reverse motion we found v"^^^ ^ lOud 
for the following parameters: amax ~ 47500, L = 8.3, 
and Ip = 160. These results suggest that the fastest 
forward swimmers will in general not move accordingly 
in the opposite direction above the transition frequency. 
Also, we note that the maximum velocities will be af- 
fected by the cross section of the filament, which in gen- 
eral will be quite smaller than the one considered here. 
Reducing the cross-section reduces the drag, and thus in- 
creases the final velocity. In particular, we note the drag 
on our model swimmer comes mostly from the tail since 
it contributes ~ 1/(1 + L^^) of the overall friction coefli- 
cient (neglecting hydrodynamic interactions). Reducing 
this quantity by a factor of 5 should yield velocities that 
are of the order of lOOx larger than its typical diffusion 
motion. Thus, it could be foreseable to use this type of 
microswimmer for perfoming tasks at the micron scale 
with high precision. 

Our results also highlight the difference for externally 
actuated microswimmers between constant driving force 



(or torque) and constant driving frequency at fixed field 
strength, since they lead to different qualitative behavior. 
In the particular case studied here, the former "ensem- 
ble" would have not led to any instability, while as shown 
here, the latter "ensemble" does exhibit a non-intuitive 
motion-reversal transition. This point is relevant because 
in most situations one deals with constant driving fre- 
quencies, and not with constant torques. Apart from 
this, our results also show that increasing the value of 
the field does not lead to large changes in the maximum 
velocity that the microswimmer can attain. Another im- 
portant point is due, and it concerns the anchoring of the 
filament. In all this work we have assumed that the semi- 
flexible polymer can relax the torsion along its contour, 
or that the anchoring is " slippy" . Although this type of 
attachment can be constructed, a more realistic scenario 
is a grafted filament completely fixed at the base. In 
this case, and assuming that the polymer cannot relax 
the stresses by rotation around its backbone, the torsion 
will not be zero and can lead to the so-called twirling- 
whirling instability [11 , and many other interesting phe- 
nomena such as plectonemes [12 . Nonetheless, one can 
estimate that this transition will occur at cutwi ~ 5w* 
for our fastest swimmer. Hence, our results should be 
clearly observed in the laboratory, even if the filament is 
not allowed to rotate freely around its backbone. 

In summary, we have introduced a simple, yet novel 
design for an externally driven microswimmer that could 
be self-assembled in the lab. The swimmer displays for- 
ward or backward motion depending on the value of the 
frequency of the actuating magnetic field. In particular 
it undergoes a transition at a critical driving frequency 
u!* above which it starts displaying stroke movements in 
the direction perpendicular to the plane of the rotating 
field. The maximum velocities achieved are well above 
the typical diffusion velocity, and thus the swimmer can 
be controlled to up to very small scales. The design pro- 
posed here is versatile and could be easily employed to 
study swarming behavior or hydrodynamic synchroniza- 
tion. Also, arrays of these swimmers could be used as 
ciliae-like biomimetic pumps in microfluidic applications. 
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